function [mse, evm_wide, evm_sub] = cal_evaluation_metrics(predict, ref)
if(var(ref)<db2pow(-100))
    fprintf('[warning][cal evm]ref power less than -100dB\n');
end
mse = mean(abs(predict-ref).^2,'all');
% EVM_{wide}=\frac{E[|\hat{y}-y|^2]}{E[y^2]}
evm_wide = mse/mean(abs(ref).^2,'all');

% EVM_{sub}=E[\frac{|\hat{y}-y|^2}{y^2}]
evm_sub = mean(abs(predict-ref).^2 ...
    ./abs(ref).^2,'all');

end